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The  concept  of  laser  cooling  (optical  refrigeration)  by 
luminescence  up-conversion  in  solids  dates  back  to  1929 
[1].  Pringsheim  recognized  that  thermal  vibrational  en¬ 
ergy  can  be  removed  by  anti-Stokes  fluorescence  if  a 
material  is  excited  with  photons  having  energy  below 
the  mean  fluorescence  energy  (see  Fig.  1).  Material  purity 
problems  prevented  observation  of  this  type  of  laser 
cooling  until  1995,  when  it  was  first  demonstrated  in 
ytterbium-doped  glass  [2].  This  was  followed  soon  after 
by  reports  of  cooling  in  dye  solutions  [3]  and  thulium- 
doped  glass  [4].  Attaining  net  cooling  in  semiconductors, 
however,  has  remained  elusive.  A  key  problem  has  been 
the  inability  of  luminescence  to  efficiently  escape  from  a 
semiconductor  due  to  total  internal  reflection  [5,6]. 
Although  the  theory  of  semiconductor  cooling  has  been 
tackled  previously  [7,8],  the  critical  issues  of  lumines¬ 
cence  trapping  and  redshifting  have  not  been  taken  into 
account.  These  processes  have  the  potential  to  frustrate 
attempts  to  achieve  semiconductor  net  cooling.  Here,  we 
resolve  this  problem  and  show  that  laser  cooling  of  semi¬ 
conductors  is  feasible. 

A  primary  advantage  of  semiconductors  compared  to 
rare-earth  doped  solids  is  their  potential  for  achieving 
temperatures  —10  K  and  below.  This  is  due  to  the  differ¬ 
ence  of  the  ground  state  populations  in  the  two  systems. 
As  the  temperature  drops  below  100  K  in  a  rare-earth 
doped  system,  the  population  at  the  top  of  the  ground 
state  manifold  dramatically  decreases,  rendering  the 
cooling  process  highly  inefficient.  This  is  a  consequence 
of  Boltzmann  statistics.  Semiconductors,  on  the  other 
hand,  obey  Fermi-Dirac  statistics,  which  keeps  the  lower 
energy  valence  band  populated  even  at  absolute  zero. 

In  this  Letter,  we  describe  laser  cooling  in  semicon¬ 
ductor  structures  allowing  for  arbitrary  external  effi¬ 
ciency.  Our  analysis,  for  the  first  time,  accounts  for  the 
luminescence  redshift  due  to  reabsorption.  We  combine 
this  theory  with  the  established  plasma  theory  of  semi¬ 
conductor  optical  absorption  including  many-body 
Coulomb  effects  and  band  filling.  The  necessary  experi¬ 
mental  conditions  for  observing  laser  cooling  in  bulk 
semiconductor  heterostructures  are  obtained. 

We  investigate  an  intrinsic  (i.e.,  undoped)  semiconduc¬ 
tor  system  uniformly  irradiated  with  a  laser  light  at 
photon  energy  hv.  We  make  the  realistic  assumption 
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that  only  a  fraction  ( r]e )  of  the  luminescence  can  escape 
the  material;  the  remaining  fraction  (1  —  r]e)  is  trapped 
and  reabsorbed,  which  causes  the  reexcitation  of  electron- 
hole  pairs  and  subsequent  reemission.  This  is  known  as 
photon  recycling.  For  steady-state  conditions  at  a  given 
temperature,  the  electron-hole  ( e-h )  carrier  density  ( N )  is 
obtained  from  the  following  equation: 

0  =  a^’N)j  -AN-  BN2  -  CN 3  +  (1  -  v  e)BN2, 
hv 

(i) 

where  1  is  the  laser  irradiance  and  a{v,  N)  is  the  interband 
absorption  coefficient.  Density  dependent  absorption  re¬ 
sults  from  (i)  Coulomb  screening  and  (ii)  band  filling.  The 
latter  effect  is  saturation  due  to  the  Pauli  exclusion  prin¬ 
ciple  (i.e.,  Pauli  band  blocking).  It  is  accounted  for  by 
taking  a{N,  hv)  =  a0(N,  hv){fv  —  fc},  where  a{)  is  the 
unsaturated  absorption  coefficient  [9].  The  bracketed 
term  is  the  band-filling  factor  defined  by  the  Fermi  dis¬ 
tribution  functions  fv  and  fc  for  the  valence  and  con¬ 
duction  bands,  respectively.  The  recombination  processes 
are  nonradiative  (AN),  radiative  (BN2),  and  Auger  (CN2). 
All  the  above  coefficients  are  temperature  dependent. 
Rearranging  the  terms  in  Eq.  (1)  yields  a  radiative  re¬ 
combination  rate  that  scales  as  rjeB.  A  similar  result  has 


FIG.  1.  Cooling  cycle  in  laser  refrigeration  of  a  semicon¬ 
ductor  in  which  a  laser  photon  with  energy  hv  is  absorbed 
followed  by  emission  of  an  up-converted  luminescence  photon 
at  hvf. 
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been  derived  by  Asbeck  [10].  It  is  important  to  emphasize 
that  rje  is  averaged  over  the  entire  luminescence  spec¬ 
trum,  i.e.,  rje  =  j  S{v)R{v)  dv /  f  R(v)dv.  Here  S{v)  is 
the  geometry-dependent  escape  probability  and  R(v)  is 
the  luminescence  spectral  density,  which  is  related  to  the 
absorption  coefficient  through  the  van  Roosbroeck- 
Shockley  relation  [11].  The  radiative  recombination  is 
obtained  from  BN 2  =  / R{v)dv  [9]. 

The  net  optical  power  deposited  in  the  semiconductor 
is  the  difference  between  the  absorbed  power  and  escaped 
luminescence: 

Pnet  =  BeBN2{hv  —  hv f)  +  ANhv  +  CN3hv  +  A  P, 

(2) 

where  we  define  an  “escaped”  mean  luminescence  pho¬ 
ton  energy  hv y  =  f  S(v)R(v)hv dv /  f  S{v)R(v)  dv.  A 
residual  heating  term  A P  accounts  for  free-carrier  ab¬ 
sorption  and  other  parasitic  absorptive  processes. 
Equation  (2)  completely  describes  laser  cooling  of  a  semi¬ 
conductor.  It  incorporates  luminescence  trapping  via  in¬ 
hibited  radiative  recombination  ( rjeB )  and  includes 
redshifting  of  the  escaped  luminescence  hPt. 

In  the  limit  of  high  external  efficiency  where  S{v)  =  1, 
Eq.  (2)  approaches  previous  results  that  assume  r]e  =  1 
and  Pf  =  Vf  [5,7,8].  Equation  (2)  indicates  laser  cooling 
occurs  when  Pnet  <  0,  which  requires  a  dominant  contri¬ 
bution  from  radiative  recombination.  Cooling  efficiency 
T]c  is  defined  as  the  ratio  of  cooling  power  density  Pc  = 
—  Pnet  to  the  absorbed  laser  power  density  ( Pabs  =  a  I  + 
A P).  Ignoring  the  A P  contributions  for  the  moment,  rjc  is 

Vc=Vi—~h  (3) 

V 

where  77 ,  =  rjeBN2/{AN  +  rjeBN2  +  CN 3)  describes  the 
“inhibited”  or  external  radiative  quantum  efficiency. 
Assuming  the  radiative  term  dominates  recombination 
(i.e.,  rjeB/C  >  N  >  A/ rjeB)  and  ignoring  band  filling, 
we  can  analytically  solve  for  optimum  carrier  density 
and  incident  laser  irradiance.  We  account  for  parasitic 
absorption  by  taking  A P  =  abI  +  afcaNI,  where  ah 
represents  background  parasitic  absorption  from  defects 
and  impurities  and  oyca  is  the  free-carrier  absorption 
cross  section.  Net  cooling  (Pnet  <  0)  occurs  if  A  <  A0, 
where  A0  defines  a  “break-even”  nonradiative  recombi¬ 
nation  rate: 


— 


ab  Y(VeB)2 
a{v)j  4  C' 


(4) 


where  we  have  defined  a  quantum  cooling  efficiency  r\q  = 
{hpf  —  hv)/hv  and  C'  =  C  +  a fcar)eB  /  a{v).  Cooling 
takes  place  for  a  laser  irradiance  in  the  range  lx<  1  < 
I2,  where  /12  =  {hprjeB/ a{v)}N2 2,  and  the  e-h  densi¬ 
ties  are 


Inclusion  of  band  filling,  however,  will  limit  the  maxi¬ 
mum  carrier  density  attainable  with  optical  pumping 
(Vmax)  that  occurs  when  fc  =  fv.  As  we  discuss  below, 
cooling  may  be  inhibited  at  low  temperature,  where  N\ 
can  approach  and  even  exceed  /Vmax .  The  parameters  B 
and  C  are  fundamental  properties  of  a  semiconductor; 
these  coefficients  can  be  calculated  and  measured  for 
bulk  or  quantum-confined  structures.  Reported  values, 
however,  vary  considerably.  In  bulk  GaAs,  for  example, 
the  numbers  range  from  2  X  10-16  <B  <1 X  10-16m3/s 
and  1  X  10~42  <  C  <  7  X  10~42  m6/s.  To  assess  the  fea¬ 
sibility  of  laser  cooling,  we  use  average  values  in  our  cal¬ 
culations:  B  =  4  X  10-16  m3/s  and  C  =  4  X  10-42  m6/s 
and  initially  ignore  background  and  free-carrier  absorp¬ 
tion.  The  break-even  nonradiative  lifetime  at  room  tem¬ 
perature  =  1/A0  is  plotted  in  Fig.  2  as  a  function  of 
r]e.  We  set  hvf  —  hv  =  kT  and  use  the  mean  lumines¬ 
cence  wavelength  of  Ay  ~  860  nm.  The  gray  area  under 
the  curve  designates  the  undesired  heating  region.  The 
horizontal  line  corresponds  to  Tnr  =  40  /rs,  which  is  the 
longest  reported  nonradiative  lifetime  in  a  GaAs/InGaP 
double  heterostructure,  although  5-10  /rs  is  typical  [12]. 
These  lifetimes  set  a  lower  limit  on  extraction  efficiency 
rje  =  15%  for  the  best  material;  in  practice  20%-30% 
may  be  required. 

Equation  (4)  suggests  increasing  the  quantum  effi¬ 
ciency  r]q  by  decreasing  the  incident  photon  energy  re¬ 
duces  the  break-even  lifetime.  At  hvf  —  hv  =  2  kT,  for 
example,  the  required  lifetime  decreases  by  a  factor  of  4. 
As  the  photon  energy  moves  into  the  Urbach  tail,  how¬ 
ever,  the  interband  absorption  drops  rapidly  and  therefore 
background  and  free-carrier  absorption  are  no  longer 
negligible.  Recently,  it  was  found  that  crfca  at  band- 
edge  wavelengths  is  much  smaller  than  previously  ex¬ 
pected  [13].  For  GaAs,  crfca  ~  10“24  m2,  which  requires 
a(v)  >  103  m_1  to  ensure  that  free-carrier  losses  are 


FIG.  2.  The  break-even  nonradiative  lifetime  as  a  function  of 
the  luminescence  extraction  efficiency  in  bulk  GaAs  at  300  K. 
This  calculated  lifetime  delineates  the  cooling  zone  and  heat¬ 
ing  zone  (shaded  area).  The  horizontal  dashed  line  corresponds 
to  Tnr  =  40  /us,  the  longest  lifetime  reported  in  a  GaAs/InGaP 
double  heterostructure. 
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negligible  compared  to  Auger  recombination  (i.e.,  C  — 
C ).  This  condition  is  satisfied  even  at  A  =  890  nm  (cor¬ 
responding  to  hvf  —  hv  =  2  kT),  where  a(v)  —  104  m”1. 
Therefore,  free-carrier  absorption  does  not  present  a  ma¬ 
jor  limitation  to  laser  cooling.  Possible  sources  for  ab  are 
(i)  impurity  absorption  in  active  and/or  cladding  layers 
and  (ii)  substrate  absorption.  It  is  also  important  that  ah 
in  Eq.  (4)  be  scaled  to  the  device  geometry:  a'b  =  ahX 
{d/L),  where  d  and  L  are  the  thicknesses  of  the  lossy  and 
active  media,  respectively.  For  an  active  medium  L  — 
1  i± m  and  substrate  thickness  of  1  cm,  this  requires 
a{ v)  >  105a/.  This,  in  turn,  demands  a'h  <  10  m”1  at 
A  =  870  nm,  and  a'h  <0.1  m”1  at  A  =  890  nm.  A  ZnS 
substrate  with  a'h  S  10“2  m”1  meets  such  a  requirement. 
A  ZnSe  substrate  having  a'h  S0.3  m”1  can  be  effective 
only  at  A  —  870  nm  [14].  We  shall  later  discuss  using  such 
substrates  to  make  index-matching  dome  lenses  to  en¬ 
hance  rje. 

Our  analysis  to  this  point  has  been  concerned  only  with 
net  cooling  from  room  temperature.  Of  key  importance  is 
determining  a  break-even  lifetime  for  cryogenic  condi¬ 
tions.  Using  the  expected  scaling  C(T)  exp(l  —  300/T) 
[15],  B  °c  T”3/2  [9,16],  keeping  rjq  —  kT/Eg,  and  ignor¬ 
ing  parasitic  losses  and  the  small  temperature  dependence 
of  the  band-gap  energy,  we  obtain 


A0(T) 

A0(300) 


/300 


300  -T\ 
T  )' 


(6) 


This  indicates  that  at  T  =  100  K,  for  example,  the  break¬ 
even  nonradiative  lifetime  is  lowered  by  —20  times  com¬ 
pared  to  room  temperature.  Furthermore,  nonradiative 
recombination  is  primarily  surface  recombination,  which 
decreases  exponentially  with  temperature  [16].  This 
makes  observation  of  laser  cooling  even  more  favorable 
at  a  lower  starting  temperature,  despite  decreased  effi¬ 
ciency  ( ~kT/Eg ).  Reduced  cooling  efficiency  at  a  lower 
temperature  is  primarily  due  to  the  smaller  probability  for 
phonon  absorption  by  electrons  and  holes.  This  is  mani¬ 
fest  in  a  temperature-dependent  exciton  linewidth  [9]: 

T(T)  =  r„  +  aT  +  yNLO(T),  (7) 


where  T0  is  due  to  an  inhomogeneous  background  distri¬ 
bution  of  impurities,  a  is  the  contribution  from  acoustic 
phonon  scattering,  y  is  the  coefficient  of  LO-phonon 
scattering,  and  NLO(T)  is  the  Bose-Einstein  phonon  dis¬ 
tribution.  Exciton-exciton  scattering  is  negligible  at  the 
densities  encountered  here  [17].  At  lattice  temperatures 
approaching  10  K,  the  acoustic  phonon  component  domi¬ 
nates  and  the  scattering  rate  becomes  comparable  to  the 
radiative  recombination  rate  (BN2).  This  means  that  cold 
exciton  recombination  occurs  before  complete  thermal- 
ization  with  the  lattice.  Hot  exciton  recombination  is  a 
similar  process  that  has  hindered  observation  of  Bose- 
Einstein  condensation  in  semiconductors.  This  problem 
has  been  significantly  alleviated  by  employing  quantum- 
confined  systems  where  relaxation  of  wave-vector  con¬ 


servation  in  the  confined  direction  reduces  a  by  nearly  3 
orders  of  magnitude  [18].  Quantum  confinement  may 
allow  realization  of  semiconductor  cooling  below  10  K. 

Band  filling  can  saturate  long  wavelength  absorption 
because  of  the  lower  density  of  states  near  the  band 
edge.  To  estimate  its  importance  in  semiconductor  cool¬ 
ing,  we  calculate  N\  using  the  same  parameters  as  in 
Fig.  2;  we  assume  an  external  efficiency  rje  =  20%  and 
temperature-independent  nonradiative  lifetime  1/A  = 
10  /rs.  The  maximum  attainable  density  is  )Vmax  —  25/Vj 
at  300  K,  dropping  to  /Vmax  —  2.5Vi  at  30  K,  with  in¬ 
sufficient  excitation  for  cooling  predicted  at  10  K.  This 
worst  case  estimate  ignores  reduced  nonradiative  recom¬ 
bination  and  overestimates  excitonic  saturation  at  low 
temperature.  A  more  accurate  model  of  band  filling  re¬ 
quires  a  rigorous  description  of  exciton  absorption  and 
luminescence  in  the  presence  of  many-body  interactions. 

We  now  analyze  the  effect  of  radiation  trapping  and 
photon  recycling  in  laser  cooling  of  GaAs.  Luminescence 
removal  is  facilitated  by  a  nearly  index-matched  dome 
attached  to  the  heterostructure  [Fig.  3(a)]  [5].  Domes 
made  from  ZnS  or  ZnSe  are  transparent  in  the  infrared 
and  present  a  minimal  heat  load  from  blackbody  radia¬ 
tion;  they  do  not  significantly  limit  the  amount  of  cooling 
that  can  be  attained.  Two  heterostructure  designs  have 
been  suggested  for  improving  the  efficiency  of  light 
emitting  diodes  as  depicted  in  Fig.  3(b).  In  Structure  I, 
both  surfaces  are  polished,  which  limits  the  radiation 
escape  cone  to  the  angle  of  total  internal  reflection.  For 
Structure  II,  the  bottom  surface  is  textured  to  statistically 
randomize  the  photon  trajectories.  More  light  rays  avoid 
the  total  internal  reflection  condition  and  escape  via 
photon  recycling.  Ignoring  Fabry-Perot  effects,  escape 
probabilities  for  both  cases  have  been  evaluated  with 
statistical  geometric  optics  using  solid-angle  averaging 
in  a  photon-gas  model  [19].  Sample  thickness  (L),  absorp¬ 
tion,  and  refractive  index  enter  the  calculation  through 


I.  Polished  U.  Textured 


FIG.  3.  (a)  A  high  index  (nd)  dome  is  bonded  to  a  semi¬ 

conductor  heterostructure  of  index  ns  consisting  of  a  GaAs 
active  layer  sandwiched  between  two  thin  layers  of  AlGaAs  or 
InGaP  for  surface  passivation,  (b)  Heterostructure  designs  with 
polished  planar  surfaces  (I)  and  a  textured  bottom  surface  (II). 
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FIG.  4.  Calculated  escape  efficiency  (top),  redshifting  of  the 
escaped  mean  luminescence  photon  energy  (middle),  and  total 
cooling  power  (bottom)  as  a  function  of  the  GaAs  thickness 
for  Structure  I  (solid  lines)  and  Structure  II  (dashed  lines) 
depicted  in  Fig.  3.  A  density  of  5  X  1017  cm-3  and  T  =  300  K 
are  assumed. 


the  previously  defined  escape  parameter  S{v).  Taking  the 
thickness  of  the  active  GaAs  as  the  independent  parame¬ 
ter,  we  analyze  these  systems  assuming  a  ZnS  dome  with 
index  of  refraction  nd  =  2.4.  ZnS  has  extremely  low 
optical  absorption  in  the  spectral  region  of  interest. 

To  evaluate  the  luminescence  spectrum,  R{v),  we  use 
an  analytical  plasma  theory  that  has  shown  good  agree¬ 
ment  with  experimental  results  for  T  >  70  K  [20].  The 
calculated  escape  efficiency  (rje)  and  shift  of  escaped 
mean  luminescence  energy  are  shown  in  Fig.  4.  Results 
are  presented  for  L  >  0.01  jjl m  where  quantum  confine¬ 
ment  of  carriers  is  not  important.  Photon  recycling  pro¬ 
vided  by  the  textured-surface  structure  results  in  much 
higher  escape  efficiency  as  the  thickness  is  reduced.  This 
enhanced  efficiency  is  obtained  at  the  expense  of  delete¬ 
rious  redshifting  of  the  luminescence  spectrum — lower 
energy  photons  suffer  less  absorption  and  can  eventually 
escape.  A  polished  (planar)  surface  exhibits  a  relatively 
flat  efficiency  up  to  thickness  1  /rm  and  the  mean  escaped 
photon  energy  remains  close  to  the  internal  value  of  hvf. 
Samples  with  thickness  L  >  0.6  jam  are  larger  than  the 
photon  mean-free  path;  both  surfaces  then  provide  essen¬ 
tially  the  same  efficiency.  A  quantity  of  interest  for  ap¬ 
plications  is  the  cooling  power.  We  define  a  normalized 
cooling  power  by  assuming  a  fixed  excitation  energy  at 
hv  =  hvf  —  kT  in  Eq.  (2): 

Pc  =  Krje(hvf  -  hv)V  «  Krje(kT  -  Ae)V,  (8) 

where  A"  is  a  normalization  constant,  A e  =  hvf  —  h Vr  is 
the  amount  of  luminescence  photon  energy  shift,  and 
V  =  Area  X  L  is  the  illumination/recombination  volume. 
We  evaluate  this  quantity  (Pc)  for  both  geometries  as  a 


function  of  the  thickness  as  shown  in  Fig.  4.  Highest 
cooling  power  and,  hence,  lowest  temperature  is  achieved 
with  a  polished  structure  of  thickness  of  1-2  jam.  If 
highest  cooling  efficiency  is  desired  and/or  a  large  rje  is 
required  to  achieve  the  break-even  condition,  a  thin  (. L  < 
500  nm)  textured-surface  structure  is  advantageous. 

In  summary,  we  have  developed  a  model  to  treat 
laser  cooling  in  semiconductors  for  arbitrary  external 
efficiency.  The  break-even  nonradiative  lifetime  and  lu¬ 
minescence  extraction  efficiency  are  obtained.  Photon 
recycling  causes  luminescence  reabsorption  and  red- 
shift,  which  may  adversely  affect  cooling  in  textured- 
surface  structures.  Highest  cooling  power  is  achieved 
with  polished  planar  structures  and  index-matching 
encapsulation. 
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